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ABSTRACT 


This  report  describes  a theoretical  invest iRation  of  the 
oscillatory  motions,  stability,  and  orientation  of  the  Moored 
Body  Scheme  A in  current  and  waves.  Linear,  first-order  equations 
for  oscillatory  motions  and  cable  tension  are  developed  for  a body 
of  revolution  moored  to  the  bottom  in  finite  depth  water  by  a 
single  restraining  cable.  Wave  exciting  forces  and  certain  hydro- 
dynamic  coefficients  are  developed  from  potential  flow  theory,  and 
non-linear  viscous  damping  effects  are  Included  in  the  formulation 
by  an  equivalent  linearization  technique.  Surge,  sway,  pitch,  and 
yaw  motion  and  tension  amplitudes  are  obtained  in  response  to 
regular  and  irregular  waves  for  a variety  of  water  depths,  cable 
lengths,  and  static  cable  tensions.  Equilibrium  orientation  in  waves 
and  currents  is  established  by  the  current  or  wave-current  inter- 
action depending  on  the  magnitude  of  body  oscillation  velocities 
relative  to  the  current  speed.  The  most  prominent  feature  of  the 
oscillatory  motion  is  the  presence  of  a pitch  resonance  occurring 
at  a frequency  of  0.33  radians  per  second. 


INTRODUCTION 


I 


1L. 


The  David  W.  Taylor  Naval  Ship  R&D  Center  (DTNSRDC)  was  requested 
to  analytically  determine  the  oscillatory  motions  and  orientation  of  a 
moored  body  under  a seaway  and  subject  to  a current.  Moored  Body  Scheme  A, 
which  Is  assumed  to  have  a circular  cross-section.  Is  statically  teth- 
ered at  a prescribed  height  above  and  parallel  to  a smooth  horizontal 
bottom  by  a single  cable  attached  to  a point  on  the  body  forward  of  the 
center  of  gravity.  The  body  may  be  subject  to  the  forces  and  moments 
induced  by  a train  of  regular,  surface  waves  propagating  through  the 
finite  depth  water  at  an  arbitrary  angle  relative  to  the  body,  the  forces 
and  moments  resulting  from  a small  steady  current,  or  the  combination  of 
both  waves  and  current. 


It  is  required  that  the  oscillatory  surge,  sway,  pitch  and  yaw 
motions,  and  the  variation  In  the  cable  tension  (roll  is  not  considered) 
be  determined  over  a range  of  sea  states  as  a function  of  water  depth, 
cable  length,  static  cable  tension,  and  wave  heading  angles.  It  is  fur- 
ther required  that  the  body  orientation  and  cable  deflection  angle  in  a 
steady  current  and  any  possible  Interaction  between  a small  current  and 
waves  be  determined. 

The  first  section  of  this  report  describes  the  theoretical  formula- 
tion and  presents  motion  results  for  the  moored  body  subject  to  incident 
waves  with  no  current  present.  The  formulation  of  the  equations  of  motion  and 
tension  fluctuations  Is  linear  and  harmonic  with  respect  to  wave  amplitude 
and  body  motions.  Heave  motion  Is  assumed  to  be  restrained  by  the  cable. 

The  equations  are  coupled  in  the  surge  and  pitch  modes  and  the  sway  and 
yaw  modes  and  contain  appropriate  linear  cable  tension  functions.  The 
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wave  exciting  forces  and  moments  and  the  added  inertia  coefficients 
which  include  the  effect  of  the  bottom  and  finite  depth  water  are  devel- 
oped from  potential  theory.  Vlacous  damping  effects,  which  are  nonlinear 
in  nature,  are  also  Incorporated  in  the  equations  by  an  equivalent  lin- 
earization technique. 

Amplitudes  of  the  surge,  sway,  pitch,  and  yaw  motion,  and  the  cable 
tension  fluctuations  in  response  to  incident  sinusoidal  waves  of  unit 
amplitude  are  presented  as  a function  of  excitation  wavelength  for 
various  water  depths,  cable  lengths,  static  cable  tensions,  and  heading 
angle  conditions.  In  addition,  the  regular  wave  responses  have  been 
used  to  obtain  the  statistical  motion  and  tension  in  an  irregular  sea 
defined  by  the  Pierson-Moskowit z Sea  Spectrum. 

In  the  second  section  of  the  report,  the  orientation  of  the  body 
in  a steady  current  is  investlca*‘ed . The  static  and  dynamic  stability 
of  the  Moored  Bodv  Scheme  A in  the  pitch  and  yaw  planes  is  determined 
based  upon  hydrodynamic  coefficients  obtained  from  model  tests  on  a con- 
figuration similar  to  Scheme  A.  Once  the  stable  pitch  and  heading  angle 
relative  to  the  current  are  known,  the  static  cable  tension  and  deflection 
angle  are  determined. 

In  the  last  section  the  Interaction  of  current  and  waves  is  briefly 
examined.  For  current  speeds  which  are  less  than  or  of  the  same  order 
as  oscillatory  velocities,  a second  order  theory  is  required  for  quanti- 
tative motion  predictions.  A qualitative  estimate  of  some  second  order 
effects  indicates  that  for  specified  current  speeds  less  than  one  knot, 
the  primary  effect  of  wave-current  interaction  is  to  align  the  average 
heading  angle  with  the  current.  When  the  wave  propagation  direction 
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relative  to  the  current  is  known,  oscillatory  motions  can  then  be  approxi- 
mated by  the  zero  current  results. 


MOTION  IN  WAVES 

EQUILIBRIUM  CONDITIONS 

The  oscillatory  motions  of  Moored  Body  Scheme  A In  response  to 
sinusoidal  surface  waves  are  first  determined  In  the  absence  of  a current. 

The  hull  Is  essentially  a body  of  revolution  with  a long  parallel  middle 
section,  hemispherical  nose,  tapered  afterbody  and  small  tail  appendages, 
as  shown  In  Figure  1.  In  the  following  analysis  the  effect  of  the 
appendages  Is  neglected.  Principal  dimensions  and  other  geometric 
characteristics  are  provided  In  Table  1 for  the  body  with  neutral  buoy- 
ancy and  60  pounds  net  buoyancy. 

Sketches  defining  principal  dimensions  of  the  Moored  Body  Scheme  A, 
the  mooring  system,  and  the  wave  train  are  shown  in  Figure  2.  The  body 
Is  moored  beneath  the  surface  In  finite  depth  water  to  a smooth,  hori- 
zontal bottom  by  a single  cable  of  negligible  mass.  The  cable  is  attached 
to  the  bottom  of  the  body  at  an  arbitrary  longitudinal  location.  A small 
net  buoyancy  keeps  the  cable  vertically  taut  and  thus  the  body  is  restrained 
at  the  prescribed  cable  length  above  the  bottom.  The  static  or  equili- 
brium cable  tension  Is  given  by 

T - B-W  (1) 

s 

where  B and  W are  the  buoyancy  and  weight  of  the  body,  respectively.  It 
Is  further  assumed  that  the  longitudinal  locations  of  the  center  of  gravity 
and  buoyancy  relative  to  the  cable  attachment  point  are  initially  established 
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Table  1 - Geometric  Character  1st Ics  of  the 
Moored  Body  Scheme  A 


Overall  length 

20 . 104 

feet 

(6 .1277m) 

1 Maximum  diameter 

1 .750 

feet 

(0.5334m) 

■ 

: 

Longitudinal  distance  from  | 

nose  to  center  of  buoyancy 

9 .340 

feet 

(2  8468m) 

1 

Height  of  center  of  gravity  J 

above  baseline 

0.6833 

feet 

(0.2083m) 

1 

1 Height  of  center  of  buoyancy  j 

1 above  baseline 

0.8750 

feet 

(0.2667m) 

1 Neutral  Buoyancy  I 

1 Mass 

86.397 

s lugs 

(1260.9kg) 

1 ’ Tngitudinal  distance  from  ] 

S nose  to  center  of  gravity 

Moment  of  inertia  about 

9 .340 

feet 

(2 .8468m) 

2 

■ 

1 

1 renter  of  gravity* 

j Moment  of  inertia  about 

2405.4 

slug  „ 
-feet"^ 

(3261  .2kg-m  ) 
2 

\ 

j y or  z axis* 

! Net  Buoyancy  60  Pounds 

2751 .0 

slug  „ 
-feet'^ 

(3729.9kg-m  ) 

i 

1 

Mass 

1 

$ 

j Longitudinal  distance  from 

84 . 53 1 

s 1 ugs 

(1233 .6kg) 

1 

nose  to  center  of  gravity 
1 Moment  of  inertia  about 

9 .383 

feet 

(2 .8599m) 

o 

1 center  of  gravity* 

2457  .3 

slug 

-feet'^ 

(333 1 . 7kg-m  ) 

N 

1 Moment  of  inertia  about  „ j 

fi  y or  z axis* 

2691 .6 

slug  „ 
-feef^ 

(3649  .3kg-m“^) 

ij 

|i  *Va  lues  for  moment  of  inertia 

distribution 

L,. 

assume 

un i f orm 

mass 

j 

J 

such  that  the  longitudinal  axis  of  the  body  la  In  equlllbrlun  parallel 
to  the  calm  water  surface  and  the  bottom.  The  zero  trim,  static  ment 


I 

^1 

J 


r 


balance  about  the  cable  attachment  point  Is 

BXg-Wx^  - 0 (2) 

where  and  are  the  longitudinal  locations  of  the  center  of  buoyancy 
and  gravity  respectively  relative  to  the  cable  attachment  point  (x^  and 
are  negative  when  the  cable  la  attached  ahead  of  the  CB  and  CG) . 

The  incident  waves  are  assumed  to  be  plane,  linear  surface  waves 
of  amplitude  A and  wavelength  propagating  at  an  angle  s relative  to 
the  equilibrium  longitudinal  axis  of  the  body  as  shown  In  Figure  2. 

In  the  finite  depth  water  the  angular  frequency  is  related  to  the  wave- 
length by^ 


CJ 


(3) 


where  g Is  the  acceleration  due  to  gravity  and  h Is  the  water  depth. 

The  wave  number  K is  defined  as  usual  by  K •=  2^/\. 

A right-handed  Cartesian  coordinate  system,  Oxyz,  is  fixed  at  the 
equilibrium  position  of  the  vehicle.  The  origin  is  located  at  a point 
directly  above  the  cable  attachment  point  on  the  body  axis.  The 
positive  X axis  Is  directed  forward  and  Is  coincident  with  the  longi- 
tudinal axis  of  the  body  at  equilibrium.  The  positive  y and  z axis 
are  directed  to  port  and  upward,  respectively.  The  polar  coordinates, 
r^ and  6,  are  defined  In  transverse  sections  along  the  length  of  the  body, 
as  shown  In  Figure  2. 
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In  the  linear  analysis  the  oscillatory  motions  of  the  body  in  response 
to  waves  are  expressed  as  small  displacements  and  rotations  from  the  equili- 
brium position  (Oxyz  system) . Body  motion  in  the  vertical  direction  is 
assumed  to  be  ideally  restricted  by  the  cable.  The  time  varying  displacements 
of  the  hull  in  the  x and  y directions  are  exprefsed  as  £^(t)  and  C,(t),  respec- 
tively, and  the  angular  rotations  about  the  y and  z axes  by  '".^(t)  and 

f;,(t),  respectively.  Roll  motion  is  not  treated  in  this  analysis, 
o 

EQUATIONS  OF  MOTION 

The  body  is  assumed  to  be  initially  at  equilibrium  in  calm  water, 
moored  at  a prescribed  height  off  the  bottom  with  zero  trim  angle.  As 
a regular  wave  train  passes  over  the  water  surface,  oscillatory  forces 
and  moments  with  the  same  frequency  as  the  Incident  wave  and  directly 
proportional  to  wave  amplitude  are  developed  on  the  body  in  one  or  all 
of  the  degrees  of  freedom  depending  on  the  particular  wave  heading 
angle.  The  body  may  respond  In  the  surge,  sway,  and  yaw  modes  to  the 
horizontal  plane  excitations  and  in  the  pitch  mode  to  the  vertical 
plane  excitations.  Heave  motion  is  restricted  by  the  cable  with  the 
oscillatory  vertical  force  balanced  by  fluctuations  In  the  cable  tension. 

As  the  body  responds  to  the  wave  excitations,  other  forces  and  moments, 
due  to  cable  tension,  inertia,  and  damping  come  into  play.  The  equations 
of  notion  are  essentially  formulated  in  a linear,  harmonic  manner, 
although  nonlinear  viscous  damping  effects  are  Incorporated  by  an 
equivalent  11  nearlzation technique . It  is  assumed  that  body  motions 
and  tension  variations  are  small  fluctuations  about  equilibrium,  hence, 

Nllne-Thomson,  "Theoretical  Hydrodynamics,"  Fifth  Edition,  Macmillan 
Company,  1968. 

★ 

Actually,  small  steady  and  higher  harmonic  wave  forces  are  also  generated 
but  are  not  considered  in  the  linear  theory. 
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only  terms  which  depend  linearly  on  wave  amplitude  and  body  motions  are 
retained  In  the  equations.  It  Is  assumed  that  no  hydrodynamic  or  Iner- 
tial forces  act  on  the  cable.  Uhen  the  wave  excitations  are  sinusoidal 


functions  of  time,  resulting  body  motions  and  tension  variations  are 
also  sinusoidal  with  the  same  frequency.  The  wave  exciting  forces, 
body  motions,  and  cable  tension  fluctuation  may  be  written,  respec- 
tively, aa 

PLCt)  = <“> 

f,(t)  =-R=.il  <‘») 

T,u;=: 


where  1"1,  2,  3,  5,  6 denotes  the  surge,  sway,  heave,  pitch,  and  yaw 
modes,  respectively,  and  , and  l,,j  are  the  complex  ampli- 

tudes of  the  wave  excitations,  body  displacements,  and  tension  fluctua- 
tions, respectively. 

The  first  order  wave  exciting  force  and  moment,  F^,  acting  on  a 
body  of  revolution  submerged  in  finite  depth  water  are  developed  In 
Appendix  A.  Three-dimensional  quantities  are  obtained  by  the  usual 
strip  approximation.  In  which  two-dimensional  sectional  forces  are 
integrated  over  the  length  of  the  body.  The  formulation  Includes  the 
effects  of  the  finite  depth  water,  the  proximity  of  the  bottom,  and 
wave  diffraction.  Explicit  expressions  for  the  components  of  the  force 
and  moment  are  given  In  Equations  (A14-a)  through  (A14-e)  of  Appendix  A. 
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The  hydrodynamic  coefficients  which  will  appear  in  the  equations 


of  motion,  are  developed  in  Appendix  B.  The  added  inertia  coefficients, 
, associated  with  body  acceleration,  are  given  In  Equations  (Bl-a) 
through  (Bl-g) . The  index,  1,  denotes  the  sense  which  the  force  or 
moment  acts,  and  the  Index,  j,  denotes  the  mode  of  motion  producing 
the  force  or  moment.  The  added  Inertia  coefficient  for  axial  motion 
Is  obtained  from  the  added  mass  of  an  equivalent  spheroid,  and  the 
coefficients  corresponding  to  transverse  body  motions  are  developed 
In  a strip-wise  fashion  from  potential  theory,  which  assumes  deep 
submergence  and  includes  the  effect  of  proximity  to  the  bottom. 

The  damping  coefficients,  , associated  with  the  body  velocity, 
are  given  in  Equations  (B3-a)  through  (B3-1)  of  Appendix  B.  The  coeffi- 
cient for  axial  motion  Is  derived  from  the  submarine  derivative,  X , 

uu 

ard  all  of  the  other  damping  coefficients  corresponding  to  transverse 
motions  are  obtained  in  a strip-wise  fashion  from  linearized  cross-flow 
drag  considerations. 

By  specifying  the  body  motions  relative  to  the  coordinate  system 

located  at  a point  directly  above  the  cable  attachment  point  on  the 

equilibrium  body  axis,  heave  is  eliminated  as  a possible  degree  of 

freedom.  To  first  order,  horizontal  and  vertical  plane  motions  are 

2 

decoupled.  The  surge  and  pitch  modes  are  coupled  , as  are  the  sway  and 
yaw  motions.  Tension  fluctuations  are  contributed  by  the  vertical 
wave  exciting  force  and  pitch  motion.  All  wave  exciting  moments  and 

‘'Newman,  J.N.,  Plala,  P.,  and  Zarnick,  E.E.,  "Pitching  Motions  of  a 
Moored  Submerged  Mine  in  Waves,"  NSRDC  Report  2151,  March  1966. 
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rotational  inertia  and  hydrodynamic  coefficients  are  computed  relative 
to  the  Oxyz  axis.  The  equations  of  motion  and  the  tension  equation  are 
then  given  as  follows. 

Surge-Pitch 

The  first  order  equations  of  motion  for  surge  and  pitch  are 

(■'v-i-AJf  -l-B,,  i +<^30%- 
( I,  + A„  ) /;  + + [ T,K,  ( I + £ ) - f ^ 

- “kK  f,  - 


♦ 

i’ 


i 


where  m is  the  mass,  I is  the  moment  of  inertia  about  the  y axis,  A,,  and 

y 11 

are,  respectively,  the  added  mass  and  damping  for  surge  notion,  and 
are,  respectively,  the  added  moment  of  Inertia  and  moment  due  to  damp- 
ing for  pitch  motion  and  the  dots  above  the  motion  variables  denote,  as 
usual,  the  time  derivative.  The  third  tern  in  the  surge  equation  is  a 
restoring  force  due  to  cable  tension.  The  third  term  in  the  pitch 
equation  and  the  coupling  terms  in  both  equations  arise  due  to  the 
attachment  of  the  cable  on  the  keel  and  the  vertical  offset  of  the 
center  of  gravity  from  the  body  axis.  It  is  assumed  that  any  hydro- 
dynamic  coupling  between  the  surge  and  pitch  modes  is  negligible. 
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Using  Equations  (4a)  and  (4b)  and  taking  appropriate  time  deriva- 
tives, the  above  equations  may  be  written  In  terms  of  the  complex 
amplitudes,  and  . 

- Fci 

f„,  + [ -<j-(  r^+A„)  -.-u.  B,, 

+ T,e.(l  5;^  = /jy 

B 


When  the  appropriate  expressions  for  the  excitation  forces  and  the 
hydrodynamic  coefficients  as  given  In  Appendices  A and  B are  substituted, 
the  above  equations  may  then  be  solved  simultaneously  for  the  complex 
amplitudes  of  surge  and  pitch  motion  per  unit  wave  amplitude.  Substitu- 
tion of  and  in  Equation  (4b)  then  provides  the  time  vary- 

ing surge  and  pitch  motion. 

Sway-Yaw 

The  first  order  equations  of  motion  for  sway  and  yaw  motion  are: 


i 


s 


A 


where  Is  the  moment  of  Inertia  In  yaw,  ®22  i’®spectively , 

the  added  mass  and  damping  for  sway  motion,  and  are,  respectively, 
the  added  moment  of  Inertia  and  moment  due  to  damping  for  yaw  motion, 

and  A,,  are  added  Inertia  coupling  coefficients  and  and  B,„  are 
zo  62  r o 26  62 

damping  coupling  coefficients.  The  third  term  In  the  sway  equation  Is 
a restoring  force  due  to  the  cable  tension.  The  choice  of  the  coordinate 
system  at  a point  above  the  cable  attachment  has  eliminated  any  cable 
restoring  moment  In  the  yaw  equation.  Coupling  Is  a result  of  both 
Inertial  and  hydrodynamic  effects  In  these  modes. 

In  terms  of  the  complex  amplitudes,  equations 

of  motion  are  written  as 


A,,)  f,,  + [ - A,  J 

- x(*J  zz 

+ A6z)  - + j^-UJ  ( 

- i-u)  1 iTot  ~ ^6 


(7) 


i 
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(8) 
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vhare  the  excitation  forces  and  hydrodynamic  coefficients  are  given  In  ) 

Appendices  A and  B. 

Cable  Tension 


In  order  that  the  body  be  completely  restrained  against  heave 
motion,  it  is  required  that  the  total  cable  tension,  which  is  the  sum 
of  the  static  and  oscillatory  components,  be  greater  than  zero  and 
less  than  a value  which  would  snap  the  cable  at  any  Instant  of  time. 

0^  Tj-t-TjCtX  Ts„„p  (9) 


If  the  total  tension  drops  to  zero,  the  cable  slackens  and  the  body  may 
move  downward.  During  a later  portion  of  the  wave  cycle,  the  body  will 
move  upward  again,  and  a large  Impulsive  force  on  the  cable  can  be 
generated  as  the  cable  becomes  taut. 

Within  the  limits  of  Equation  (9) , the  vertical  force  must  be 
balanced  by  both  the  static  and  time  varying  components  of  the  cable 
tension.  Statically,  the  balance  Is  provided  by  Equation  (1).  The 
oscillatory  component  Is  given  by 

and  in  terms  of  the  complex  amplitudes 


(10) 


where  A^j  and  are  the  added  Inertia  and  damping  providing  a vertical 
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fore*  due  to  pitch  motion.  Equation  (9)  rcqulrea  that 


I \ c,-  ) 


, whichever  is  smaller. 


NONLINEAR  DAMPING  AND  THE  SOLUTION  OF  THE  EQUATIONS 

The  restoring  forces  and  moments  affecting  the  moored  body  can 
Introduce  resonant  motion  when  the  damping  is  small.  As  is  typical 
for  a linear  system,  the  amount  of  damping  determines  the  magnitude 
of  the  motion  near  the  resonant  frequency.  The  primary  source  of  damping 
for  a deeply  submerged,  oscillating  body  is  the  nonlinear  viscous  drag, 
and  although  this  force  may  be  small,  it  can  be  of  significance  if  a 
motion  resonance  occurs  In  the  frequency  range  of  excitation.  In  Appen- 
dix B a method  is  shown  whereby  the  nonlinear  viscous  drag  is  approximated 
by  a pseudo-linear  function  of  velocity  which  may  be  incorporated  in  the 
linear  equations  of  motion.  The  damping  coefficients,  • temain 
amplitude  dependent,  however,  and  the  equations  of  motion  are  solved  by 
an  Iterative  procedure. 

When  the  equations  of  motion  are  strictly  linear,  that  is,  all  of 
the  coefficients  are  functions  of  frequency  only,  solutions  are  ordin- 
arily obtained  in  the  form  of  response  amplitude  operators,  which  are 
Che  motion  amplitudes  per  unit  wave  amplitude,  I 'c  lA  . The 
response  amplitude  operators  are  Independent  of  wave  amplitude,  and 
the  absolute  amplitude  of  motion,  ^ fo.  | , Is  obtained  by  multiplying 

the  response  by  an  arbitrary  wave  amplitude.  In  the  case  where  the 
damping  coefficients  are  amplitude  dependent,  solutions  may  be  obtained 
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in  the  form  of  a response  amplitude  operator,  except  that 
Is  no  longer  Independent  of  wave  amplitude. 

From  Appendix  B,  the  equivalent  pseudo-linear  form  for  the 
damping  coefficients  is  given  by 


By  = 


where  the  constants  depending  upon  geometry  and  mode  of 

motion,  as  shown  In  Equations  (B3-a)  through  (B3-1) . In  terms  of  the 
notion  per  unit  wave  amplitude,  | ^ 1 /A  , the  damping  coefficient 

is 


In  the  Iterative  method  used  to  solve  the  equations  of  motion, 
initial  values  for  the  damping  coefficients  B^^  are  determined  from 
Initial  arbitrary  choices  for  /^y//A  with  a given  value  for  the 
wave  amplitude,  A.  These  coefficients  are  then  substituted  Into  the 
equations  of  motion,  which  are  solved  for  new  values  of  / I ^ ^ 

The  new  motion  responses  are  then  used  to  redefine  the  damping  coef- 
ficients, and  the  process  is  repeated  until  a convergent  solution  is 
obtained.  Since  the  damping  coefficient  is  a function  of  wave  amplitude, 
a distinct  response,  , is  determined  for  each  specified  wave 

amplitude. 

The  damping  term  appearing  in  the  equations  of  motion  is  also 
proportional  to  Cp,  the  crossflow  drag  coefficient  for  a circular 
cylinder.  This  coefficient  is  difficult  to  estimate  since  it  is  to 
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some  extent  an  unknown  function  of  Reynolds  number,  wave  frequency,  and 
amplitude.  The  coefficient  may  vary  from  0.4  to  1.2  depending  on  whether 
the  flow  la  completely  turbulent  or  laminar. 

In  the  following  theoretical  results,  the  product,  C^A,  has  been 
retained  as  an  arbitrary  constant  i^hlch  determines  the  magnitude  of  the 
damping  coefficient.  It  la  understood  that  and  A may  each  be  varied 
Independently  within  appropriate  limits. 

An  example  of  the  resonant  motion  and  the  effect  of  nonlinear 
damping  on  the  Moored  Body  Scheme  A is  shown  In  Figure  3.  The  pitch 
amplitude  per  unit  wave  amplitude  In  degrees  per  foot  is  plotted  versus 
excitation  wavelength.  The  body  Is  moored  by  a 20-foot  cable  at  a point 
2 feet  ahead  of  the  center  of  gravity.  The  water  depth  is  120  feet,  the 
wave  heading  angle  is  180  degrees,  and  the  static  tension  is  60  pounds. 

A pitch  resonance  occurs  at  a wavelength  of  about  1100  feet,  which 
according  to  Equation  (3)  corresponds  to  an  angular  frequency  of  0.325 
radlans/sec.  The  family  of  curves  at  the  resonance  corresponds  to  differ- 
ent values  for  the  product,  C^A.  Away  from  resonance,  the  curves  in  the 
family  merge  into  a single  curve,  thus  Indicating  the  significance  of 
the  damping  only  about  resonance. 

From  Figure  3,  the  peak  magnitude  of  resonance  may  be  plotted  against 
C,^A.  As  shovm  In  Figure  4 the  response  at  resonance  is  roughly 
proportional  to  the  inverse  square  root  of  the  damping  factor,  C^A.  The 
absolute  value  of  pitch  motion  has  then  been  plotted  against  wave 
amplitude  for  different  choices  of  C^,  and  this  is  shown  in  Figure  5. 

Other  modes  of  motion  would  be  similarly  affected  about  resonance 
by  a damping  variation,  however,  the  surge,  away,  and  yaw  resonances 
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Cq  • DRAG  COEFFICIENT  FOR  A CIRCULAR  CYLINDER 
A - WAVE  AMPLITUDE  IFEET) 


WAVELENGTH  (FEET! 


Fijiurc  i - Effect  of  Non-linear  Dampinu  on  Pitch  Response  in  Regular  Waves 


Figure  4 Effect  of  Damping  on  Figure  5 Effect  of  Wave  \nipli 

the  Peak  Magnitude  of  the  Pitch  tude  and  Drag  (’<K*fficient  on 

Response  at  Resonance  in  Regular  Pitch  Angle  at  Resonance  in 

Waves  Regular  Waves 


occur  ac  longer  wavelengths,  as  will  be  shown  in  the  next  section. 


As  a final  coement  on  the  nonlinear  damping,  the  damping  coeffi- 
cients developed  In  Appendix  B are  formulated  from  cross-flow  drag 
considerations,  and  the  velocities  In  the  drag  terms  should  Ideally 
reflect  the  relative  velocity  between  the  body  and  fluid.  In  the 
analysis,  however,  the  wave  orbital  velocity  has  been  neglected. 


MOTION  IN  REGULAR  WAVES 

The  motion  amplitudes  of  the  Moored  Body  Scheme  A in  response  to 
regular  waves  have  been  obtained  from  the  theory  developed  in  the 
previous  section  for  a range  of  water  depths,  cable  lengths,  static 
tensions  and  wave  heading  angles.  Principal  dimensions  of  Scheme  A 
are  provided  In  Table  1.  The  mooring  cable  is  attached  to  the  keel 
of  the  vehicle,  two  feet  ahead  of  the  center  of  gravity  and  all  body 
motions  are  determined  relative  to  the  coordinate  system,  Oxyz , fixed 
on  the  equilibrium  body  axis  directly  above  the  cable  attachment  point. 

It  was  specified  that  Scheme  A would  always  be  moored  at  least 
100  feet  below  the  calm  water  surface  and  at  least  20  feet  above  the 
bottom.  Water  depth  and  cable  length  are  permitted  to  vary  independently 
within  these  limits.  A nominal  static  cable  tension  of  60  pounds  and 
minimum  cable  length  of  20  feet  were  specified.  The  effect  on  motion 
and  tension  fluctuations  for  variations  from  these  nominal  values  Is 
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also  examined. 
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The  presence  of  restoring  forces  and  moments  In  the  equations  of 
motion  means  that  resonances  will  be  Introduced  In  the  motion  responses. 
The  natural  frequencies  for  the  various  modes  are  given  by 
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The  natural  frequencies  are  tabulated  in  Table  2 for  different  values 
of  the  static  cable  tension,  cable  length,  and  vertical  center  of  gravity. 

Coupling  between  the  modes  also  means  that  a resonance  in  one  mode 
may  be  evident  in  another  node.  This  is  particularly  true  of  the  sway- 
yaw  modes.  Although  the  yaw  equation  contains  no  restoring  moment,  the 
strong  coupling  due  to  Inertial  and  hydrodynamic  effects  means  that  the 
sway  resonance  is  significantly  reflected  in  the  yaw  response.  The 
pitch  resonance  is  also  reflected  to  some  degree  in  the  tension  response. 
Coupling  between  the  surge  and  pitch  nodes  is  relatively  weak  although 
the  surge  resonance  can  affect  the  pitch  response  if  the  cable  length 
becomes  very  small. 

At  frequencies  less  that  about  0.25  radians/sec  resonances  are  of  little 
consequence  in  a real  seaway  condition,  since  there  is  little  energy 
below  this  frequency  even  in  severe  sea  states.  From  Table  2 for  a 
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Table  2 - Resonant  Frequencies  of  the  Moored  Body 
Scheme  A 


Stat ic 

Cable 

Tension 

Cable 

Lenpth 

Vert  ica 1 
Center  of 
Gravity 

Resonant  Frequencies 
in  radians/second 

in  feet 

in  feet 

in  feet 

Pitch 

Sur  pe 

Sway-Yaw 

60 

10 

-0.1917 

0.326 

0.266 

0.187 

•• 

20 

• t 

0 .325 

0.188 

0.132 

50 

ft 

0.325 

0.119 

0.084 

• ( 

CO 

" 

0.325 

0.0 

0.0 

120 

10 

-0 .1917 

0.340 

0.381 

0.266 

•• 

20 

f t 

0.339 

0.269 

0.188 

•• 

50 

" 

0.338 

0.170 

0.119 

♦ ♦ 

ao 

tt 

0 .338 

0.0 

0.0 

20 

10 

-0.1917 

0.317 
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nominal  static  tension  of  60  pounds  and  cable  length  greater  than  20 
feet,  pitch  Is  the  only  mode  where  resonant  motion  Is  significantly 
excited  by  a real  seaway. 

Figures  6 through  10  show  the  effect  on  the  motion  and  tension 
responses  of  water  depth  variation  with  constant  depth  of  submergence. 

The  amplitude  of  motion  or  tension  per  unit  wave  amplitude  Is  plotted 
as  a function  of  excitation  wavelength.  The  water  depth  Is  Increased 
from  120  to  400  feet  and  the  cable  length  Is  correspondingly  increased 
from  20  to  300  feet  such  that  a 100-foot  submergence  depth  is  maintained. 
The  minimum  submergence  depth  (100  feet)  regardless  of  water  depth 
provides  a worst  case  condition  for  all  quantities.  As  Indicated  in 
the  figures,  the  120-foot  water  depth  is  a worst  case  for  surge,  sway, 
and  yaw  motion,  and  the  400-foot  depth  (essentially  infinite  depth) 

Is  a worst  case  for  pitch  and  tension  fluctuations. 

The  responses  were  obtained  with  a nominal  static  cable  tension  of 
60  pounds.  Surge  and  pitch  motions  are  maximum  In  head  waves  (B“180 
degrees) ; sway  and  yaw  motions  are  maximum  in  beam  waves  (P**90  degrees) , 
so  results  are  provided  only  for  these  heading  angles.  Tension  fluctua- 
tions are  mostly  unaffected  by  heading  angle. 

Except  for  surge,  which  Is  given  for  the  zero  damping  case,  the 
responses  for  the  other  modes  of  motion  were  obtained  for«the  damping 
factor,  CjjA"1.0.  With  the  exception  of  pitch,  all  resonances  occur  at 
very  long  wavelengths  (low  frequencies),  and  hence  the  damping  variation 
Is  not  of  great  significance  In  the  wavelength  range  Indicated  In  the 
figures.  The  effect  of  damping  variation  for  pitch  motion  Is  shown  In 


Figure  3. 
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Figure  8 - Effect  of  Water  Depth  on  Sway  Motion  in  Regular  Waves 


Figure  9 - Effect  of  Water  Depth  on  Yaw  Motion 
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From  Table  2 the  pitch  resonant  frequency  remains  essentially 
unchanged  as  cable  length  Increases.  Due  to  the  increased  water  depth, 
however,  the  resonance  is  shifted  toward  longer  wavelengths,  as  noted 
in  Figure  7.  In  fact,  the  peak  for  the  200  and  400-foot  depths  are 
shifted  beyond  the  wavelength  range  shown.  The  cable  length  increase 
does  shift  the  surge,  away,  and  yaw  resonances  to  lower  frequencies, 
and  these  peaks  are  well  beyond  the  wavelength  range  shown  in  Figures  6, 
8,  and  9.  The  tension  responses  are  shown  in  Figure  10.  Tension  is 
primarily  determined  by  the  vertical  wave  excitation  force,  although 
some  coupling  from  pitch  motion  is  evident  at  the  pitch  resonant  wave- 
length. 

Figures  11  through  15  show  the  effect  of  cable  length  variation  at 
a constant  water  depth  of  120  feet  and  nominal  static  cable  tension  of 
60  pounds.  All  responses,  except  surge,  were  obtained  with  the  damping 
factor,  CpA*=1.0.  Generally,  as  the  cable  length  Increases,  pitch 
motion  and  tension  fluctuations  Increase,  and  horizontal  plane  motions 
decrease  with  the  resonance  being  shifted  to  lower  frequencies. 

Figures  16  through  19  show  the  effects  of  static  cable  tension 
variation  for  constant  water  depth  of  120  feet  and  cable  length  of  20 
feet.  As  before,  all  responses,  except  surge,  were  obtained  with  a 
damping  factor  Cj^A-l.O.  An  Increase  in  the  static  tension  shifts  all 
resonances  towards  higher  frequencies  and  thus  Increases  the  motion 
amplitudes  at  frequencies  above  resonance.  In  the  first  order  theory, 
tension  fluctuations  are  relatively  unaffected  by  the  static  tension 
and,  hence,  are  not  shown. 
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Fijiure  1 1 - Effect  of  Cable  Length  on  Surge  Motion  in  Regular  Waves 


h = 120  FEET 
li  = 180  DEGREES 


hg  = 50  FEET 


Tj  = 60  POUNDS 


hg  = 10' 


hg  .20 


600  800 
WAVELENGTH  (FEET) 


Figure  I 2 — Effect  of  Cable  Length  on  Pitch  Motion  in  Regular  Waves 
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Figure  1 3 — Effect  of  Cable  Length  on  Sway  Motion  in  Regular  Waves 
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Figure  14  — Effect  of  Cable  Length  on  Yaw  Motion  in  Regular  Waves 
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Figure  1 8 - Effect  of  Static  Cable  Tension  on  Sway  Motion  in  Regular  Waves 


Figure  19  - Effect  of  Static  Cable  Tension  on  Yaw  Motion  in  Regular  Waves 
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The  motion  responses  shovm  In  Figures  6 through  19  are  independent 
of  wave  amplitude,  except  around  the  resonances,  and  the  pitch  motion 
is  the  only  mode  in  which  the  resonant  frequency  la  really  high  enough 
to  fall  well  within  the  range  of  realistic  wave  excitation.  It  Is 
required,  however , that  the  absolute  amplitudes  of  the  motion  and  hence, 
wave  amplitude,  are  small  enough  such  that  the  basic  linear  assumptions 
of  the  formulation  are  satisfied.  For  example,  surge  or  sway  amplitudes 
which  are  greater  than  10  feet  when  the  cable  length  is  only  20  feet 
would  be  approaching  unacceptable  limits  of  the  linear  theory. 

Roughly,  the  linear  theory  is  expected  to  provide  reliable  results 
when  the  amplitude  of  the  cable  deflection  angle  from  the  vertical  and 
the  pitch  and  yaw  angle  amplitudes  are  less  than  30  degrees. 

As  mentioned  previously,  a condition  of  the  linear  formulation  is 
chat  the  tension  fluctuation  be  small  compared  to  the  static  cable  tension. 
With  a static  tension  of  60  pounds,  this  condition  is  clearly  marginal. 

As  an  example,  consider  the  amplitude  of  tension  fluctuations  as  given  by 
Figure  10  for  the  400-foot  water  depth  and  wavelength  of  about  600  feet. 

The  tension  fluctuation  response  is  about  17  pounds/foot,  and  for  a wave 
of  even  a few  feet  amplitude  the  absolute  amplitude  of  the  tension 
fluctuation  is  of  the  order  of  the  static  tension.  Furthermore,  the 
wave  amplitude  does  not  have  to  be  very  large  before  the  vertical 
restraint  condition  of  Equation  (11)  is  violated. 

The  static  cable  tension  could  be  increased  to  relieve  this  situation, 
however,  this  would  also  have  the  adverse  effects  of  shifting  all  of  the 
motion  resonances  towards  higher  frequencies,  as  Indicated  in  Table  2. 
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The  fact  that  tension  fluctuations  may  not  necessarily  be  small 
compared  to  the  static  tension  for  the  conditions  specified  in  this 
analysis  is  a limitation  of  the  linear  analysis.  The  large  tension 


fluctuations  can  Introduce  harmonic  components  into  the  motion  responses, 
however,  a consistent  second  order  theory  would  be  required  to  quanti- 
tatively predict  the  effect. 

I 
I 

MOTION  IN  A SEAWAY 

A statistical  representation  of  the  motions  of  Scheme  A in  irregular 

seas  may  be  determined  in  a manner  following  the  procedure  of  St.  Denis 

3 4 

and  Pierson  and  summarized  by  Frank  and  Salvesen  • Briefly,  with 

the  assumptions  of  stationary,  Gaussian  processes  for  the  wave  excita- 
tions and  motion  responses,  the  energy  associated  with  the  response  of 
the  body  in  random  waves  may  be  represented  by  the  linear  superposition 
of  responses  to  a ntimber  of  sinusoidal  waves. 

The  significant  amplitude  of  motion  in  a particular  mode  is  given  by 


fo 


(12) 


where  B^ui)  is  the  amplitude  of  the  response  in  a given  mode  to  a sinusoidal 
wave  of  unit  amplitude,  and  S(u,)  is  the  energy  in  the  sea  which  is  associ- 
ated within  an  infinitesimal  frequency  band  at  the  frequency,  u.  In  this 
analysis,  it  has  been  assumed  that  the  seaway  is  unidirectional  and 

fully  developed,  and  can  therefore  be  statistically  described  by  the  one 

3 

parameter,  Plerson-Moskowitz  Spectrum  of  the  form 
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2 2 

where  A^“0.0081g  , B^«33.561/hj^^2  , and  significant  wave- 

height  in  feet.  A plot  of  the  Plerson-Moskowltz  Spectrum  (Equation  (13)) 

is  shown  in  Figure  20  for  several  values  of  Other  characteristics 

4 

of  the  Spectrum  are  given  in  Table  3 . The  significant  wavehelght  is 
the  average  of  the  highest  one-third  of  the  waves  and  is  a measure  of 
the  sea  state.  Subject  to  the  condition  of  linear  waves,  it  is  assumed 
that  the  Pierson-Moskowltz  Spectrum,  as  given  in  Equation  (13)  , is 
independent  of  water  depth,  although  the  relationship  between  frequency 
and  wavelength  given  by  Equation  (3)  is  still  applicable. 

The  method  just  outlined  to  determine  the  statistical  motion  in 
irregular  waves  is  applicable  for  strictly  linear  responses.  For  a 
given  response  amplitude  operator,  which  is  Independent  of  wave  ampli- 
tude, and  a specified  sea  spectrum,  the  significant  amplitude  of  motion 
is  uniquely  defined  for  a given  significant  wavehelght.  In  the  present 
analysis  the  nonlinear  damping  and  pseudo-linearization  technique  used 
to  solve  the  equations  of  motion  result  in  motion  responses  which  are 
amplitude  dependent  near  resonances.  However,  even  though  the  response 
amplitude  operators  are  not  uniquely  defined  independent  of  wave  amplitude, 

3 

St.  Denis,  M.  and  W.S.  Pierson,  "On  the  Motion  of  Ships  in  Confused 
Seas,"  Trans.  SHAME,  Vol.  61,  1953. 

4 

Frank,  W.  and  N.  Salvesen,  "The  Frank  Close-Fit  Ship-Motion  Computer 
Program,"  NSRDC  Report  No.  3289,  June  1970. 


Table  3 - The  Pierson-Moskowi tz  Sea  Spectrum 
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we  may  use  the  statistical  method  to  define  an  envelope  for  the  signi- 
ficant motion  and  tension  amplitudes.  As  shown  In  Figure  3,  the  product 
CjjA  provides  limits  on  the  statistical  amplitudes  at  each  significant 
wavehelght . 

The  response  amplitude  operators  like  those  shown  In  Figures  6 
through  19,  together  with  the  Pierson-Moskowitz  Sea  Spectrum,  as  defined 
in  Equation  (13),  have  been  used  In  Equation  (12)  to  obtain  the  signifi- 
cant amplitude  of  the  motions  and  tension  fluctuation  as  a function  of 
significant  wavehelght  (sea  state).  Results  are  shown  in  Figures  21 
through  32  for  the  various  water  depth,  cable  length,  and  heading  angle 
conditions  previously  considered.  In  all  cases,  the  static  cable  ten- 
sion was  60  pounds.  With  the  exception  of  surge,  results  were  obtained 
with  damping  factor  limits  of  0.2_^  C^A  < 2.0.  The  envelope  defined  by 
these  limits  Is  Indicated  by  the  shaded  region  in  the  figures.  Surge 
results  are  shown  only  for  the  zero  damping  case,  since  in  practice 
damping  In  this  mode  is  very  small.  Response  amplitude  operators  were 
computed  between  frequencies  of  0.25  and  1.4  radlans/sec.  As  shown  in 
Figure  20,  for  significant  wavehelghts  less  than  30  feet,  virtually 
all  of  the  wave  energy  la  contained  in  this  range. 

Figure  21  shows  the  significant  amplitude  of  pitch  obtained  with 
response  amplitude  operators  defined  for  the  various  values  of  the  damp- 
ing coefficients,  Cj^A.  For  sea  states  with  significant  wavehelght  less 
than  about  17  feet  the  spectral  density  at  low  frequencies  (long  wave- 
lengths) Is  insignificant  and,  hence,  resonant  motion  is  not  excited.  In 
higher  sea  states  the  spectral  density  at  lower  frequencies  is  no  longer 
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negligible,  and  the  long  waves  tend  to  generate  motion  about  the  resonance 
where  the  response  is  amplitude  dependent.  Reasonable  limits  for  the  damp- 
ing factor  may  then  be  regarded  as  defining  an  envelope  for  the  statistical 
motion  amplitudes.  Figures  23  and  28  show  the  pitch  envelope  for  water 
depth  and  cable  length  variations  within  the  range,  0.2  < ♦ 2.0. 

Figures  24,  25,  29,  and  30  show  the  significant  amplitudes  for 
sway  and  yaw.  Envelopes  for  these  modes  are  much  narrower  than  the 
pitch,  since  the  resonances  occur  at  very  low  frequencies  and  within 
this  range  of  excitation  there  is  very  little  change  in  the  response 
amplitude  operators  due  to  the  damping  variation. 

In  Figures  22  and  27  are  shown  the  significant  amplitude  of  surge 
motion.  These  results  indicate  no  envelope  since  only  the  zero  damping 
response  amplitude  operator  was  used. 

Figures  26  and  31  show  the  significant  amplitude  of  the  tension 
fluctuations.  These  responses  are  relatively  unchanged  by  the  damping 
variation,  which  would  appear  only  through  pitch  coupling. 

Figure  32  shows  the  effect  of  heading  angle  on  the  various  motion 
and  tension  significant  amplitudes.  A significant  wavehelght  of  25 
feet  was  chosen  and  the  shaded  region  Indicates  the  damping  factor 
limits  0.2^CpA  ^2.0. 

At  high  sea  states  the  significant  amplitudes  of  motion  shown  in 
Figures  21-32  indicate  magnitudes  which  exceed  the  valid  range  of  the 
linear  theory.  In  particular  a condition  for  the  linear  formulation 
of  the  equations  of  motion  and  restraint  of  the  body  in  the  heave  mode 
is  that  the  cable  tension  fluctuations  not  exceed  the  static  cable 
tension  (Equation  (11)).  Figure  26  shows  that  in  a seaway  the  signflcant 
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Figure  21  - Effect  of  Damping  on  Pitch  Motion  in  Irregular  Waves 
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Figure  23  — Effect  of  Water  Depth  on  Pitch  Motion  in  Irregular  Waves 
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Figure  29  — Effect  of  Cable  Length  on  Sway  Motion  in  Irregular  Waves 
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Figure  30  - Effect  of  Cable  Length  on  Yaw  Motion  in  Irregular  Waves 
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SIGNIFICANT  AMPLITUDE  TENSION 
(POUNDS) 


SIGNIFICANT  AMPLITUDE 


amplitude  of  the  tension  fluctuations  can  exceed  the  nominal  60  pound 


static  cable  tension  at  even  moderate  sea  states  (h  > 10  feet  for 
the  400  foot  cable  length).  The  motion  responses  shown  in  Figures  21-32 
have  been  extended  into  the  high  sea  state  range  in  excess  of  the  tension 
limits  shown  in  Figure  26  as  though  the  linear  formulation  was  valid. 

In  the  high  sea  states,  however,  the  curves  should  be  regarded  as  only 
qualitatively  indicative  of  the  vehicle  statistical  motions. 

BODY  ATTITUDE  IN  A STEADY  CURRENT 

The  orientation  of  the  Moored  Body  Scheme  A in  a steady  current 
is  examined  in  this  section.  Possible  equilibrium  pitch  and  heading 
angles  are  determined  from  static  nonlinear  pitch  and  yaw  moment  equa- 
tions, respectively.  To  obtain  the  stable  equilibrium  body  orientations 
in  the  current,  the  dynamic  stability  of  the  body  in  the  pitch  and  yaw 
modes  at  each  static  equilibrium  position  is  developed  using  the  stability 
equations  derived  in  Appendix  C.  The  static  and  dynamic  equations  for 
vertical  and  horizontal  plane  motion  utilize  hydrodynamic  coefficients 
based  upon  model  tests  performed  on  a configuration  similar  to  the  Moored 
Body  Scheme  A. 

STATIC  EQUILIBRIUM 

A definition  sketch  of  the  moored  body  in  a current  is  shown  in 
Figure  33.  The  body  has  a prescribed  net  buoyancy  and  is  assumed  to  be 
trimmed  at  zero  pitch  in  no  current.  It  is  assumed  that  the  cable  is 
massless  and  that  no  hydrodynamic  forces  act  on  it.  A coordinate  system 
Ox'y'z'  is  fixed  on  the  body  axis  at  a point  directly  above  the  cable 
attachment  point,  as  shown.  Orientation  in  the  vertical  and  horizontal 
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planes  will  be  treated  Independently.  When  the  body  Is  pitched  at  an 
angle  relative  to  the  current,  Z4  * ccSjc.  and 

denote  the  negative  of  the  current  velocity  components  In  the  x'  and  z* 
directions,  respectively.  When  the  body  Is  yawed  at  an  angle 

/ 

relative  to  the  current,  = UeCCC.|^^  and  denote  the 

negative  of  the  current  velocity  components  In  the  x'  and  y'  directions, 
respectively. 

Relative  to  the  body  fixed  system  located  above  the  cable  attachment 
point,  the  static  pitch  moment  as  a function  of  and  the  static  yaw 
moment  as  a function  of  are,  respectively, 

AV^C'^o)  = A\  K)  + 'Z  (.<o) 

Vr(A)=A'(A)f 

where  AIK)  and  are  the  static  hydrodynamic  pitch  and  yaw 

moments  respectively  about  the  center  of  gravity,  and  and 

^ (S<r)are  the  static  hydrodynamic  forces  In  the  z' and  y*  directions, 
respectively.  The  static  moment  due  to  the  cable  tension  and  center 
of  gravity  are  Included  In  the  pitch  moment  equation. 

The  static  hydrodynamic  forces  and  moments  developed  on  the  Moored 
Body  Scheme  A have  been  estimated  based  upon  experimental  data  for  a 
shortened  version  of  the  body.  The  data  was  modified  to  account  for 
the  Increased  length  of  the  Scheme  A configuration,  and  the  resulting 
estimated  static  coefficients  are  given  In  Table  4.  The  notation  follows 
that  commonly  used  in  standard  submarine  equations.^  The  coefficients 

^Certler,  M.  and  Hagen,  G.,  "Standard  Equations  of  Motion  for  Submarine 
Simulation,"  DTNSRDC  Report  2510,  June  1967. 


(14a) 

(14b) 
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Table  4 - Non-Dimensional  Static  HyUrodymamic  ^ 

Coefficients  for  the  Moored  Body  Scheme  A 


^ , Y'^  ^ -0.00676 

, -N'^  0.00571 

-=  -0.05771 

M^^  , -N'^  , -0.00680 

♦Moments  referred  to  the  center  of  gravity.  Equilibrium 
angle  of  attack  is  zero . 


Table  5 - Non-Dimensional  Stability  Derivatives 
for  the  Moored  Body  Scheme  A* 


z; . 

Y'v  . 

-0.01090 

m;  , 

1 

II 

0.00435 

-0.00899 

< = 

-Ni  = 

-0.00081 

Z'  . 
q 

-0.00377 

M'q  = 

N'r  = 

-0.00181 

z'i,  - 

-Vf  - 

-0.00068 

“i  = 

nJ.  = 

-0.00082 

I'y  , - 0.00086 

m'  ^ m/  , 0.01075 

x'j,  ^ X(j/L  . -0.09948 


* 

Moments  referred  to  0 system.  Equilibrium 
angle  of  attack  is  zero. 
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are  expressed  In  nondlmenslonal  form  and  are  given  relative  to  the  body 
center  of  gravity.  Because  of  the  geometric  symmetry  of  the  body, 
vertical  plane  coefficients  are  equally  applicable  (with  appropriate 
sign  changes)  for  the  horizontal  plane. 

The  static  hydrodynamic  force  and  moment  as  a function  of  angle  of 
attack  were  originally  measured  over  a + 15  degree  range.  The  resulting 
data  were  then  fitted  with  a second  order  polynomial  In  slno^o  and  cosu^^, 
the  coefficients  of  which  are  the  static  coefficients  given  In  Table  4, 
and  modified  for  the  Scheme  A configuration. 

The  static  pitch  and  yaw  moment  relative  to  the  point  on  the  body 
axis  directly  above  the  cable  attachment  as  given  by  Equations  (14a)  and 
(14b),  respectively,  may  then  be  rewritten  In  nondlmenslonal  form  In  terms 
of  the  static  coefficients  given  In  Table  4, 

_ ( aVu,  OCS  I SmxJ.  ) 

" , '-T;  ^ -*^330.1  (15a) 

at  *- 

^ CCSi9,.5trrV$c  ■¥  5m  Pt  1 Swv|3ol  ) 

^ ^ f ! r I 

- ^ (V^  '♦"'lylvl  I Swr\  Pol  ) (15b) 

In  Figure  34,  the  nondlmenslonal  static  pitch  moment  given  by 

Equation  (15a)  Is  plotted  as  a function  of  angle  of  attack,  for  several 

values  of  the  current  speed.  The  hydrodynamic  coefficients  are  given 

in  Table  4,  T "60  pounds,  h„  " <*>  , and  the  vertical  center  of  gravity 
s a 

is  located  a nominal  0.19  feet  below  the  body  axis. 


L 
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Due  to  the  symmetry  of  the  vertical  and  horizontal  plane  hydro- 
dynamic  coefficients  Equations  (15a)  and  (15b)  become  Identical  at  large 
current  speed,  so  that  the  Infinite  speed,  static  pitch  moment  curve 
also  represents  the  static  yaw  moment,  which  Is  Independent  of  speed. 

Possible  equilibrium  pitch  and  heading  orientations  relative  to  the 
current  may  occur,  respectively,  when  the  total  static  pitch  and  yaw 
moment  are  zero.  Referring  to  Figure  34  In  the  vertical  plane,  for 
current  speeds  less  than  two  knots,  the  total  moment  Is  zero  only  at 
zero  angle  of  attack  (pitch  angle).  For  current  speeds  greater  than 
approximately  two  knots,  three  equilibrium  pitch  angles  are  possible. 

At  Infinite  current  speed  0 and  + 23  degree  pitch  angles  are  potential 
equilibrium  orientations.  In  the  horizontal  plane,  heading  angles  of 
0 and  + 23  degrees  relative  to  the  oncoming  current  are  possible 
equilibrium  orientations  at  all  current  speeds. 

DYNAMIC  STABILITY 

In  order  to  determine  whether  the  static  equilibrium  conditions 

obtained  in  the  previous  section  are  stable  body  orientations  in  a 

current,  it  is  necessary  to  determine  the  dynamic  stability  of  the 

moored  body  at  each  possible  static  equilibrium  pitch  and  heading 
orientations.  In  Appendix  C linear  equations  for  dynamic  stability 

of  the  Moored  Body  Scheme  A In  a current  are  developed.  The  motion 

equations  are  formulated  for  small  perturbations  In  the  vertical  and 

horizontal  planes  about  arbitrary  equilibrium  pitch  and  heading  angles, 

respectively,  and  the  dynamic  stability  in  each  plane  Is  determined 

from  the  roots  of  the  governing  characteristic  equation. 

The  hydrodynamic  force  and  moment  coefficients  used  In  the  deter- 
mination of  the  stability  criteria  are  the  stability  derivatives,  which 
represent  the  slopes  of  the  force  or  moment  with  respect  to  a given 
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motion  variable  at  a prescribed  angle  of  attack  (or  drift  angle  p,.  ). 

The  stability  derl\ratlves  at  zero  angle  of  attack  for  a configuration 
of  the  Moored  Body  vhlch  Is  shorter  In  length  are  a’^allable  from 
a previously  conducted  program  of  experirents. 

These  coefficients  were  modified  to  account  for  the  Increased  length 
of  the  Scheme  A configuration,  and  then  adjusted  for  the  fact  that  the 
origin  of  the  body  fixed  coordinate  system  was  located  directly  above 
the  cable  attachment  point.  These  derivatives  are  provided  in  Table  5. 

Since  the  derivatives  are  taken  about  zero  angle  of  attack,  a quantita- 
tive description  of  the  dynamic  stability  is  possible  only  at  ^ and 

p - O • However,  based  upon  estimated  values  for  the  derivatives  at 
other  angles,  some  qualitative  predictions  can  be  made  regarding  the  dynamic 
stability  at  equilibrium  pitch  and  heading  angles  other  than  zero. 


Vertical  Plane 


In  Appendix  C the  basic  equations  for  coupled  surge  and  pitch  motion 
of  the  moored  body  in  a current  are  given  by  Equations  (C3)  and  (C4)  . This 
system  of  equations  results  in  a fourth  order  characteristic  equation, 
the  roots  of  which  determine  the  nature  of  the  coupled  surge  and  pitch 
motion. 


A simpler  condition  for  vertical  plane  stability  is  obtained  from 
the  equation  for  uncoupled  pitch  motion.  A condition  for  dynamic  stability 
in  the  pitch  mode  which  neglects  the  effect  of  surge  coupling  but  includes 


I 


(non-rlgorously)  the  restoring  moment  due  to  the  vertical  center  of  gravity 
is  given  by  Equation  (C7)  of  Appendix  C, 
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For  an  equilibrium  pitch  angle  of  zero  (a^  - 0),  with  - 60  pounds, 

Zq  “ -0.19  feet,  and  given  from  Table  5 vertical  plane  stability  is 

obtained  when  U <2.4  knots, 
o 

Figure  34  Indicates  that  for  current  speeds  roughly  greater  than 

two  to  three  knots,  a static  equilibrium  occurs  at  angles  of  attack  other 

than  zero;  specifically  + 23  degrees  at  infinite  current  speeds.  The 

stability  derivative  M ' is  only  available  at  a “0,  its  value  at  other 

w o 

equilibrium  angles  may  be  qualitatively  estimated  from  Figure  34,  however. 

The  slope  at  a =0  of  the  nondimenslonal  total  moment  curve  for  infinite 
o 

current  speed  given  in  Figure  34  represents  the  stability  derivative  M 

w 

At  a =0  this  slope  is  equal  to  the  value  of  M ' given  in  Table  5,  and  is 

o V 

a positive  quantity.  For  static  angles  of  attack  greater  than  approximately 

10  degrees,  the  slope  of  the  infinite  speed  total  moment  curve  becomes 

negative,  therefore  M ' is  negative  and  from  Equation  (16)  dynamic  stability 

w 

is  obtained  whenever  is  negative. 

The  vertical  plane  stability  characteristics  for  the  Moored  Body 
Scheme  A with  T » 60  pounds,  z = -0.19  feet  and  h ■=  OO  may  be  summarized 

8 0 O 

as  follows 

(a)  For  U < 2.4  ^<-nots  the  body  is  dynamically  stable  at  zero  equilibrium 

o 

pitch  angle. 

(b)  For  U >2.4  knots  the  body  becomes  dynamically  unstable  at  zero 

o 

pitch  angle  but  becomes  dynamically  stable  at  plus  and  minus  equilibrium 
pitch  angles  which  deviate  from  zero  angle  with  increasing  current  speed. 

At  infinite  current  speed  the  body  is  stable  at  + 23  degree  pitch  angles. 
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Horizontal  Plane 


From  Figure  34  the  static  yaw  moment  acting  on  the  moored  body  is 
given  by  the  curve  for  infinite  speed  which  is  zero  at  heading  angles  of 
+ 23  degrees  and  zero  degrees  relative  to  the  oncoming  current.  In  order 
to  determine  the  orientation  the  body  actually  assumes,  it  is  necessary 
to  determine  the  dynamic  stability  at  these  possible  equilibrium  headings. 
The  expressions  for  coupled  sway  and  yaw  motion  of  the  moored  body 
in  a current  are  given  by  Equations  (C-10)  and  (C-11)  of  Appendix  C.  This 
system  of  equations  results  in  a cubic  characteristics  equation  (Equa- 
tion C-12)  the  roots  of  which  determine  the  nature  of  the  coupled  sway 
and  yaw  motion.  Using  the  horizontal  plane  stability  derivatives  for 
zero  drift  angle  (6  “ 0)  as  given  in  Table  5,  Equation  (C-12)  becomes 


JL 


(17) 


where  o'  “ s is  the  nondlmenslonal  root,  and  s = d/dt. 
o 

For  a nominal  static  cable  tension  of  60  pounds  and  minimum  cable 

length  of  20  feet,  the  stability  requirements  given  by  Equation  (C-13) 

of  Appendix  C are  satisfied  when  the  current  speed  ‘ 1.75  knots,  Wlien 

the  cable  length  is  Increased  to  50  feet,  the  body  is  dynamically  stable 

at  P =0  for  U < 1.10  knots.  For  infinite  cable  length  Equation  (17) 
o o 

reduces  to  a quadratic,  the  third  teim  of  which  is  negative,  thus  indi- 
cating the  body  to  be  unstable  at  all  current  speeds. 

Specific  values  for  the  stability  derivatives  are  not  available  for 
the  + 23  degree  equilibrium  headings.  However,  the  critical  derivative 
affecting  the  stability  characteristics  is  N^'  and  this  coefficient  may 
be  estimated,  as  in  the  case  of  vertical  plane  stability,  from  the  slope 


of  the  Infinite  speed  moment  curve  of  Figure  34.  From  Figure  34  it  is 
noted  that  the  slope  of  the  total  moment  at  “ + 23  degrees  is  of 
opposite  sign  to  the  slope  at  6^  ■ 0.  Substitution  of  a large  positive 
value  for  N^'  into  Equation  (C-12)  means  that  the  third  term  of  Equation  (17) 
becomes  positive  for  all  combinations  of  cable  tension  and  length  and 
current  speed.  Anticipated  changes  in  the  other  stability  derivatives 
at  “ 23  degrees  are  not  expected  to  significantly  alter  the  charac- 
teristic equation.  With  all  of  the  terras  in  Equation  (17)  of  the  same 
sign  and  Equation  (C-13)  satisfied,  the  moored  boy  is  stable  at  6^  “ + 23 
degrees  for  all  conditions. 

The  horizontal  plane  stability  characteristics  of  the  Moored  Body 
Scheme  A may  be  summarized  as  follows: 

(a)  The  body  is  dynamically  stable  at  heading  angles  of  + 23 
degrees  relative  to  the  oncoming  current  for  all  current  speed, 
tension,  and  cable  length  conditions. 

(b)  The  body  is  also  dynamically  stable  at  a zero  degree  heading 
for  current  speeds  less  than  1.75  knots  when  the  cable  length  is 
minimum  (20  feet).  As  the  cable  length  increases  the  stable 
speed  range  at  zero  heading  is  reduced  until  the  body  is  dynami- 
cally unstable  at  infinite  cable  length  for  all  current  speeds. 

(c)  Although  the  moored  body  with  finite  cable  length  is  stable 
at  low  current  speeds  at  zero  heading  angle,  stability  at  the  + 23 
degree  headings  is  stronger,  and  it  is  expected  the  body  would 
assume  these  heading  orientations  subject  to  large  perturbations. 
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CABLE  ORIENTATION 


An  analysis  of  the  dynamic  stability  of  the  Moored  Body  Scheme  A 
in  a current  indicates  that  for  low  current  speeds  the  body  assumes  a 
zero  pitch  attitude  and  possible  + 23  degrees  heading  angles  relative  to 
the  oncoming  current.  With  this  given  orientation  of  the  body  relative 
to  the  current,  the  static  cable  deflection  angle  and  cable  tension  may 
be  determined. 

The  total  lateral  force  as  a function  of  static  heading  (drift) 
angle  is  given  by 

^ ~ ^9  Si/fv^o  ” ^IV)  Si/TV I ^‘''’^^•1  ') 

where  the  coefficients  are  given  in  Table  4. 

A force  balance,  as  shown  in  Figure  33,  provides  expressions  for 
the  static  cable  deflection  and  static  cable  tension. 

Ts  = [(B-w)‘  + K(  p./  ] 

With  6^  ” 23  degrees,  Figures  35  and  36  show  the  static  cable  deflection 
and  tension,  respectively,  for  several  values  of  the  buoyancy  minus  weight. 
For  current  speeds,  less  than  about  one  knot,  the  cable  deflection  Is  less 
than  20  degrees  and  static  tension  Is  virtually  the  same  as  the  net  buoyancy. 


|- 
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WAVE-CURRENT  INTERACTIONS 


f 


In  the  first  part  of  this  report,  the  oscillatory  motions  of  the 
Moored  Body  Scheme  A In  response  to  exciting  waves  were  determined  in 
the  absence  of  a current.  In  the  second  section  the  stability  and 
orientation  of  the  body  In  a steady  current  In  the  absence  of  wave 
excitations  have  been  determined.  Clearly,  these  two  cases  must  repre- 
sent the  extremes  when  the  simultaneous  effect  of  current  and  waves  are 
present.  We  now  examine  the  behavior  of  the  body  when  it  Is  subject  to 
waves  and  a steady  current. 

In  order  to  quantitatively  determine  the  Interaction  of  .raves 
and  current,  when  current  speed  and  oscillatory  velocities  are  of  the 
same  order,  a formulation  of  the  equations  of  motion  accurate  to  second 
order  Is  necessary.  This  Is  beyond  the  scope  of  the  present  analysis. 
However,  a brief,  and  by  no  means  rigorous,  examination  of  certain 
second  order  effects  has  provided  some  qualitative  predictions  of  the 
wave-current  coupling. 

Briefly,  the  second  order  equations  In  the  yaw  plane  were  formulated. 
The  form  Is  essentially  that  of  Equations  (C8)  and  (C9)  in  Appendix  C 

except  that  the  nonlinear  force  and  moment,  Y ■ ivlvl  and  N i ivlv! 

V I v I ' ' V v 1 ' ' 

are  included,  and  where  the  lateral  velocity  v appears,  the  relative 
velocity  ( U(,  j'-rv.  ) would  be  substituted.  Also,  the  wave 

exciting  forces  appear  on  the  right-hand  side  of  the  equation.  Steady 
second  order  forces  and  moments,  due  to  the  wave  excitation,  have  been 
neglected  since  the  body  Is  submerged  at  least  100  feet  below  the  surface. 

' By  taking  the  Fourier  components  of  each  term  in  the  second  order 
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equations,  separate  equations  for  the  steady  forces  and  the  first  order 
harmonic  forces  are  obtained.  Solution  of  the  steady  equations  pro- 
vides the  mean  orientation,  and  solution  of  the  first  harmonic  equation 
provides  the  oscillatory  motions.  Results  of  this  analysis  Indicate 


(a)  The  steady  yaw  moment  equation  reduces  to  Equation  (ISb)  when 


oscillation  motions  are  zero,  and  the  first  harmonic  equations  reduce 
to  Equations  (7)  and  (8)  when  current  speed  Is  zero. 

(b)  The  combination  of  steady  current  and  oscillatory  body  motion 
in  the  cross-flow  drag  terms  provides  an  additional  steady  force  and 
moment  which  tends  to  orient  the  body  Into  the  current.  When  current 
speed  Is  much  larger  than  oscillatory  velocities,  this  steady  force  is 
proportional  to  the  square  of  body  motion  amplitude.  When  the  oscillatory 
velocity  is  much  larger  than  current  speed,  the  steady  force  Is  propor- 
tional to  the  product  of  current  speed  and  oscillatory  notion  amplitude. 

(c)  The  forces  and  moments  due  to  lift  Vy  V and  V , respec- 
tively, provide  linear  damping  In  the  first  harmonic  equations  propor- 
tional to  the  product  of  current  speed  and  oscillatory  velocity.  This 
is  the  predominant  damping  tern  when  the  current  speed  Is  large  compared 
to  oscillatory  velocities.  When  oscillatory  velocities  are  large,  the 
nonlinear  damping  proportional  to  the  oscillation  velocity  squared  is 
predominant . 

Basea  upon  results  of  the  second  order  analysis,  the  behavior  of 


the  Moored  Body  Scheme  A in  response  to  various  combinations  of  waves 
and  current  may  be  summarized  as  follows: 


1.  When  the  vehicle  ie  subject  to  a current  with  no  wave  excitation. 
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stable  heading  angles  of  + 23  degrees  relative  to  the  current  direction 
are  usually  attained,  although  under  certain  conditions  of  current 
speed,  cable  length  and  tension,  a zero  degree  heading  can  also  be 
achieved.  When  the  current  speed  Is  less  than  about  2.4  knots,  the 
vehicle  Is  stable  at  a zero  degree  pitch  angle. 

2.  When  very  small  wave  excitations  are  superimposed  on  the 
current,  the  body  assumes  mean  displacement,  heading,  and  pitch  angles 
determined  by  the  steady  current.  The  oscillatory  motions  may  be 
regarded  as  small  perturbations  about  the  equilibrium  position. 

The  equations  governing  oscillatory  motion  contain  linear  damping, 
the  coefficient  of  which  Is  directly  proportional  to  the  current 
speed. 

3.  As  oscillatory  velocity  Increases,  but  remains  less  than  current 
speed,  steady  forces  and  moments,  resulting  from  the  nonlinear  drag, 
tend  to  reduce  the  stable  heading  angles  relative  to  the  current 
direction.  Both  the  linear  damping  due  to  lift  and  the  nonlinear 
damping  due  to  cross-flow  drag  appear  In  the  equations  for  oscillatory 
motion . 

4.  When  current  speed  Is  of  the  same  order  or  greater  than 
oscillatory  velocities,  the  vehicle  assumes  a heading  angle  into  the 
current,  and  the  nonlinear  drag  is  the  predominant  damping  effect  in 
the  equations  of  oscillatory  motion.  The  oscillatory  equations  are 
approximated  by  Equations  (5)  through  (8). 

5.  When  the  body  Is  excited  by  waves  with  no  current,  there  are 
no  steady  forces  generated  and  the  vehicle  assumes  no  preferred 


direction. 
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When  the  body  Is  simultaneously  subject  to  current  and  waves, 
depending  upon  current  speed,  wave  amplitude,  and  the  angle  between 
current  and  wave  propagation  direction,  the  vehicle  may  assume  any 
one  of  the  above  conditions.  However,  for  specified  current  speeds 
less  than  1 knot«  condition  4 Is  probably  the  most  likely  situation. 
From  the  oscillatory  motion  results  of  the  first  section,  the  motion 

„ s « ♦ 

velocities  ^ least  of 

order  ujA.  With  constant  wave  slope  of  A > 1 knot 

and  f f k r .and  jL  ' 

In  summary,  when  the  body  Is  simultaneously  subject  to  waves 
and  a small  current,  the  oscillatory  motions  are  represented  by  the 
zero  current  condition.  The  only  effect  of  current  then  Is  to  orient 
the  body  so  that  the  mean  heading  and  pitch  angles  are  directed  into 
the  current. 
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CONCLUSIONS 

Based  upon  the  results  attained  from  a theoretical  Investigation 
of  the  motions  of  the  Moored  Body  Scheme  A In  waves  and  current , the 
following  conclusions  are  drawn: 

1.  In  the  absence  of  any  wave  disturbance,  the  body  Is  dynamically 
stable  at  a zero-degree  pitch  angle  for  current  speeds  less  than  2.4  knots. 
For  current  speeds  greater  than  2.4  knots  the  stable  pitch  angle  shifts 
away  from  zero  and  approaches  + 23  degrees  as  the  current  speed  becomes 
extremely  large.  Tlie  vehicle  Is  stable  at  + 23  degrees  heading  angles  rela- 
tive to  the  Incident  current  at  all  current  speeds,  and  a stable  zero  degree 
heading  angle  may  also  be  achieved  for  limited  combinations  of  current 


speed,  cable  length  and  static  cable  tension.  For  specified  current 
speeds  less  than  1 knot,  the  static  cable  deflection  from  the  vertical 
is  less  than  10  degrees,  and  the  static  cable  tension  is  not  slgnfl- 
cantly  changed  from  the  prescribed  net  buoyancy  of  60  pounds. 

2.  When  wave  Induced  oscillatory  velocities  are  small  compared 
to  current  speed,  oscillatory  motions  take  place  about  a mean  orienta- 
tion defined  by  the  steady  current. 

3.  For  current  speeds  less  than  1 knot,  oscillation  velocities 
are  typically  greater  than  the  current  speed.  In  this  case  oscillatory 
motions  are  approximated  by  the  zero  current  formulation,  the  only 
effect  of  wave  coupling  being  to  align  the  vehicle  so  that  the  mean 
pitch  and  yaw  angles  are  Into  the  current. 

4.  The  zero  current  formulation  of  wave  induced  oscillatory 
motions  indicates  that  the  mooring  system  Introduces  resonances  In  all 
of  the  modes  of  motion.  Pitch  motion,  however,  is  the  only  mode  where 
the  resonant  frequency  is  high  enough  to  be  grossly  excited  in  a real 
seaway.  The  pitch  resonance  occurs  at  approximately  0.33  radlans/sec. 

Due  to  the  nonlinear  character  of  the  damping  forces,  a precise  pre- 
diction of  the  motion  near  resonance  is  not  possible.  However,  an 
envelope  for  the  motion  defined  by  reasonable  limits  of  the  damping 
coefficient  indicates  pitch  motion  amplitudes  in  excess  of  30  degrees 
in  sea  states  6 and  7. 

5.  A minimum  submergence  depth  of  100  feet  and  infinite  water 
depth  provide  a worst  case  for  pitch  motion  and  cable  tension  fluctua- 
tions. The  100-foot  depth  of  submergence  and  minimum  120-foot  water 

\ 

i 

I 
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depth  provide  a worst  case  for  horizontal  plane  motions.  An  Increase 
in  static  cable  tension  tends  to  Increase  all  motion  quantities  by 
shifting  resonances  to  higher  frequencies. 

6.  For  a nominal  static  cable  tension  of  60  pounds,  the  signi- 
ficant amplitudes  of  the  motions  and  tension  for  the  worst  case  condi- 
tions are: 

Sea  State  $ 5 C^-/3<  13 

Surge  and  sway  do  not  exceed  two  feet.  Pitch  and  yaw  angles 
do  not  exceed  10  and  1 degrees,  respectively.  Cable  tension 
fluctuations  can  exceed  static  cable  tension. 

Sea  States  6 _ and  7 

The  cable  tension  fluctuations  exceed  the  static  cable  tension 
(the  significant  amplitude  of  is  about  170  pounds  at  ^ 

25  feet)  and  the  body  can  experience  motion  in  the  heave  mode. 

The  linear  formulation  of  the  equations  of  motion  are  no  longer 
strictly  valid  in  this  range.  Qualitatively  all  motion  increases 
rapidly  with  increasing  sea  state,  and  the  pitch  motion  may  be  in 
excess  of  30  degrees  at  * 20  feet. 
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APPENDIX  A 


WAVE  EXCITING  FORCES 


In  the  main  body  of  this  report  the  linear  equations  of  motion 
in  five  degrees  of  freedom  are  developed  for  a slender  body  of  revolution 
tehered  near  a horizontal  bottom  and  subject  to  the  forces  Induced  by 
incident  plane  surface  waves.  On  the  right-hand  side  of  the  motion 
equations  appear  the  complex  amplitudes  of  the  wave  exciting  forces, 

, where  i ■ 1,2, 3, 5, 6 correspond  to  the  surge  sway,  heave,  pitch 
and  yaw  modes,  respectively. 

In  this  section  expressions  for  the  first  order  wave  exciting 

forces,  which  Include  body  diffraction  and  bottom  effects,  are  developed 

from  finite  depth,  two-dimensional,  potential  flow  theory.  Three- 

dimensional  forces  are  obtained  in  a strlpwlse  sense  by  Integrating 

transverse  sectional  forces  over  the  length  of  the  body. 

We  assume  that  the  fluid  is  inviscid,  incompressible,  and  Irrota- 

tlonal,  the  incident  waves  are  linear  and  the  wavelength  large  compared 

to  the  transverse  dimensions  of  the  body  (2^R  /X « < 1),  and  the  bodv  is 

o 

deeply  submerged  such  that  it  behaves  as  a weak  scatterer  of  the  incident 

* 

waves . 

If  p denotes  the  pressure  on  the  body  surface,  S , the  force  and 

o 

moment  acting  on  the  body  can  be  expressed  by^ 


(Al-b) 


This  means  that  the  disturbance  created  by  wave  diffraction  does  not 
alter  the  free  surface. 


r»v**>*» 


where  n is  a unit  normal  vector  on  and  positive  into  the  body  surface, 

and  T is  a vector  from  the  origin  of  the  coordinate  system  to  a point 
o 

on  the  body  surface. 

Since  the  fluid  is  assumad  Irrotatlonal , a velocity  potential, 
i'(x,  y,  z,  t),  the  gradient  of  which  represents  the  velocity  field, 
in  the  fluid,  can  be  Introduced.  The  pressure  in  the  fluid  is  then 
related  to  the  velocity  potential  by  Bernoulli's  Equation  for  unsteady 


flow,  which  is  given  as 


-il 


where  the  static  head  and  velocity  squared  terms  have  been  neglected. 

Tor  a body  in  waves  and  restrained  from  motion,  the  velocity  poten- 
tial consists  of  two  parts;  one  representing  the  disturbance  generated 
by  the  undisturbed  incident  wave,  ♦j,  and  the  other  due  to  the  scattering 
or  diffraction  of  the  incident  wave  from  the  body, 

+ (A3) 

The  velocity  potential  for  plane,  progressive,  linear  waves  propagating 
In  a finite  depth  fluid  can  be  expressed  by^ 


rm  ». re  relationship  between  the  wavelength  and  frequency  is  given  by 
1)  In  the  main  text. 


rigid  body  surface  the  normel  velocity  component  of  the  fluid 
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^ ^ ior.)  = 0 


thus  on  S 


Inn  =-ll 


|is  = - ^ - ;..t  J 


yj^  in,  cosp -t- la^stnp -f- Dji’nk  KCj+k)j 


where  n2»  and  n^  are  the  components  of  n along  the  x , y , and  z 
directions,  respectively.  Since  n^^ « n2  or  n^  for  a slender  body  except 
near  the  body  ends,  we  can  neglect  the  contribution  of  n^^  in  the  above 
expression  and  obtain 


= -J^cgsAK(jft'c) 

GO  cosk  K h J J 

X ^ xn.j,  <3i/rv ^ -f  n^tanii  K k^)  j (A5) 

Formally,  is  obtained  by  the  solution  of  a boundary  value  problem, 
where  in  addition  to  the  usual  free  surface,  radiation,  and  bottom  condi- 
tions, must  also  satisfy  Equation  (A5)  on  the  body  surface,  S^.  Lee 
and  Newman^  have  shown  for  the  infinite  depth  case  that  an  approximate 
expression  for  "Jp  can  be  obtained  in  terms  of  the  two-dimensional  potentials 
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for  forced  harmonic  motion,  4'3(y*z)»  where 


= n and  ^ n, 


on  S 


Applying  this  same  approach  to  the  finite  depth  case,  4>q  can  be 
given  to  leading  order  in  Ky  and  Kz  by 

J„=  +f.s] 


(A6) 


where 


y = 


_ cos 


k K^c 


cosh.  Kk 


and 


^ _ si/r^l<  Kc 
C05K.  K ^ 


The  velocity  potentials,  <}i2(y»z)  and  4i3(y,z),  correspond  to  the 
two-dimensional  flow  disturbances  generated  by  the  forced  oscillation 
of  the  body  in  the  horizontal  and  vertical  directions,  respectively. 
When  the  body  Is  deeply  submerged  such  that  the  generated  disturbances 
do  not  disturb  the  free  surface,  Lee  and  Newman^  showed  that  ^2  and  i(>3 
can  be  expressed  In  the  form 

( r,  0 ; O = - ^£1^  [ -f  ) A-  ] 

M 71'  f'  ' 

^ (.r.e;  O - - f a,,uV(--] 


(A7) 

(A8) 


Lee,  C.M.,  and  J.N.  Newman,  "The  Vertical  Mean  Force  and  Moment  of 
Submerged  Bodies  Under  Waves,"  Journal  of  Ship  Research,  Vol.  15, 
No.  3,  September  1971. 


72 


where  S^(x)  la  the  sectional  area  and  822 (*)  and  a^^Cx)  are  the  sectional 
added  masses  for  motion  In  the  sway  and  heave  directions,  respectively. 

The  added  masses  for  a circular  cylinder  translating  In  the  trans- 
verse directions  near  a rigid  wall,  Is 


f I + fip-  ] 1 


(A9) 


Using  the  expressions  for  a22t  a^^t  '^'2  ^3»  approximate  expression 


for  ❖p  Is  given  by 


(AlO) 


I.  ) (-  ^V  un/5  CC^  _ j ] 

When  the  body  is  a slender  body  of  revolution,  the  surface  integral 
can  be  written  In  terms  of  the  cylindrical  coordinates  as 


n 


jj.:.  ,h  = p.f-  /?C.) 


The  normal  vector  can  also  be  written 


X. 


(All) 


and 


K = (n,n^,n,)  ^ -Si^iO  ) 

Cx»r  = ( C 'k  •.  ' ) 


(Al2-a) 


(Al2-b) 


kennard , E.  H.,  "Irrotatlonal  Flow  of  Frictionless  Fluids,  Mostly  of 
Invariable  Density,"  NSPDC  Report  2299,  1967. 


SubsClCuClon  of  Equations  (A2) , (A3) , and  (All)  Into  Equations 
(Al-a)  and  (Al-b)  provide  the  components  of  the  wave  exciting  force 


and  moment  as  follows: 


/•n 


h -- 


f>  (cix.k\x)  (do  a: 


(A13) 


where  1 > 1,  2,  and  3 denotes  forces  In  the  x,  y,  and  z directions  and 

1 “ 5 and  6 denotes  the  pitch  and  yaw  moments,  respectively. 

The  expressions  for  and  given  by  Equations  (AA)  and  (AlO) , 
respectively,  and  the  surface  normal  components  given  by  Equations 
(A12-a)  and  (A12-b)  may  then  be  substituted  In  Equation  (A13) . With 

y “ r cos  6 and  z ■ r sin  6,  the  Integrand  can  be  expanded  in  powers  of 

KR(x)  on  the  body  surface,  and  the  integration  around  the  circular  cross- 
section  of  leading  order  terms  may  be  accomplished  expllcltely.  Resulting 
expressions  then  Involve  lengthwise  Integrations  which  may  be  computed 
once  the  radius  is  specified  along  the  length  of  the  body.  The  expressions 
for  the  wave  exciting  forces  and  moments  in  five  degrees  of  freedom  are: 


r3 


(A14-a) 


i 


(AlA-b) 


(A14-C) 


= M S e [ 1 1 S,(x)  X fi  + 1^]  e ' ' 

= p>cl^y.^^  i‘^[  U s,w f 


(A14-d) 


where  the  real  part  of  all  the  above  expressions  is  to  be  taken,  and 
and  6 arc  defined  by  Equations  (A6) . 
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appendlx  b 

HYDRODYNA-MIC  coefficients 
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APPENDIX  B 


HYDRODYNAMIC  COEFFICIENTS 

The  hydrodynamic  coefficients  appearing  in  the  equations  of  motion 
are  the  added  inertia  coefficients,  , associated  with  the  acceleration 
of  the  body  and  the  damping  coefficients,  j * associated  with  the  velocity. 
Except  for  the  axial  coefficients,  Aj^^^  and  other  coefficients  are 

obtained  from  a strip  approximation  in  which  sectional  quantities  are 
computed  and  then  integrated  over  the  length  of  the  body  to  provide 
three-dimensional  results. 

ADDED  INERTIA  COEFFICIENTS 

The  added  inertia  coefficient  or  added  mass  for  axial  motion,  A^^^, 
may  be  estimated  by  assuming  the  vehicle  to  be  an  equivalent  spheroid. 

g 

Hence,  the  added  mass  is  given  by 


where 


- / 


- VA'' 


All  other  added  inertia  terms  can  be  obtained  from  the  strip  approxi- 
mations as  follows: 


Lamb,  H.  Hydrodynamics,  6th  ed.,  Cambridge  1932, 
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^21  ~ j ^22^^  ) 

A 

(Bl-b) 

Ass  = ) 
A 

^33  = ] 

A 

(Bl-c) 

a 

(Bl-d) 

(Bl-e) 


(Bl-f) 


(Bl-g) 


where  integration  is  carried  over  the  length  of  the  body  and  a22(x)  and 
a^jCx)  are  the  added  masses  of  a section  at  location  x in  the  sway  and 
heave  modes,  respectively.  If  we  assume  that  the  sections  are  deeply  sub- 
merged such  that  any  forced  motion  does  not  disturb  the  calm  water  surface, 
then  the  added  masses  for  circular  cylinders  moving  near  a rigid  wall  are 
given  by^ 


f-  2. -I 


(B2) 


where  R(x)  is  the  circular  cylinder  radius  at  location  x. 


DAMPING  COEFFICIENTS 


When  the  body  is  submerged  far  enough  below  the  water  surface  such 
that  no  disturbance  on  the  calm  water  surface  is  generated  by  forced  body 
motion,  the  so-called  wavemaking  damping  or  potential  forces  depending 
linearly  on  body  velocity  are  zero.  There  are,  however,  viscous  forces 
acting  on  the  body  which  typically  depend  upon  the  square  of  the  velocity, 
and  therefore  cannot  be  substituted  directly  into  the  linear  equations  of 
motion. 


1 
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9 

We  can,  however,  use  the  method  of  equivalent  linearization  to 
replace  the  nonlinear  drag  by  pseudo-linear  terma  such  that  the  equations 
of  motion  may  be  solved. 

The  axial  viscous  drag  can  be  written  in  terms  of  the  submarine 

derivative  X given  in  Reference  5 as 
uu  “ 


If  we  attempt  to  cast  the  drag  in  the  formi  ic  is  clear  that  the 

coefficient  is  a time  varying  function  of  velocity.  The  equivalent 

q 

linearization  method  replaces  f’V  the  time  averaged  quantity. 


3it 


5oi 


, such  that  the  damping  coefficient  may  be  written  as 


(B3-a) 


where  |Coil  is  the  amplitude  of  the  displacement  and  the  factor,  8/3”,  is 
the  coefficient  of  the  first  term  in  the  Fourier  expansion  of  . 

The  damping  coefficient,  is  now  time  independent,  although  still 

dependent  upon  the  amplitude  of  the  motion.  The  coefficient  may,  however, 
be  substituted  into  the  equation  of  motion,  with  an  arbitrarily  chosen 
amplitude,  and  the  equation  then  Iterated  until  a convergent  solution  for 
f,0i  is  obtained. 

The  other  damping  coefficients  may  be  derived  from  simple  cross-flow 
drag  consldera'‘lons  and  the  strip  approximation.  The  crosr-  flow  drag  per 
unit  length  on  a circular  section  has  the  form 


Fp/L  - 1/2  pD(x)  v^ 


9 

Thomson,  W.  T. , "Mechanical  Vibrations,"  Prentice-Hall,  Inc.,  Englewood 
Cliffs,  N.J.,  1963. 


where  the  Is  the  transverse  velocity  of  the  section  of  diameter, 

• • • • 

D(x),  and  may  be  associated  with  ^3»  xCs.  or  xCg.  The  symbol,  C^, 

Is  the  drag  coefficient  of  a circular  cylinder.  Casting  the  drag  In  terms 
of  the  coefficients,  B , replacing  | C. | by  8m  ^ , and  Integrating 

over  the  length  of  the  body,  the  damping  coefficients,  t expressed 
as 


(B3-b) 

^33 

(B3-c) 

1 t,  k\  \ ) 

A 

\,= 

J 

(B3-d) 

/ j 

];tU  kjlXlcix 
"l 

L 

:^TT  ' 

•'L 

(B3-e) 

^^<.6  " ^ jT  Co  1 ^00  i 
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APPENDIX  C 


r 


EQUATIONS  FOR  DYNAMIC  STABILITY  OF  THE  MOORED  BODY  SCHEME  A 
In  the  following  section,  the  criteria  for  horizontal  and  vertical 
plane  dynamic  stability  of  the  moored  body  In  a steady  current  are 
developed  from  the  linearized  equations  of  motion.  The  equations  are 
formulated  relative  to  a body  fixed  coordinate  system  where  the  origin 
Is  located  on  the  body  axis  above  the  cable  attachment  point.  The  form- 
ulation and  notations  used  follow  the  standard  equations  for  submarine 
motion  given  In  Reference  5 with  the  exception  that  the  positive  z axis 
Is  tc.ken  to  be  upward  and  the  positive  y axis  to  port  as  shown  in  Figure  33. 

tt 

The  equations  are  written  for  small  perturbations  from  the  equili- 
brium position.  It  Is  assumed  that  at  equilibrium  the  cable  is  vertical 
but  the  body  may  have  an  equilibrium  pitch  and  heading  angle  relative  to 
the  current  denoted  by  and  6^,  respectively,  as  shown  In  Figure  33. 

The  quantities  u,  v,  and  w are  the  perturbation  velocities  along  the  body 


VERTICAL  PLANE 


In  the  vertical  plane  the  axial  and  pitch  equations  are  coupled, 


w and  are  kinematically  related  by 


Axial  Force: 


Is  3^c  'i/t 


Pitch  Moment: 


Ij  ^ f il  - m j.£,  t,  u'  - ^ = pi  i- 

+ tp.  -t-  Pi  Liu/  (, 

+ -"yjc- J ‘Is  + i 


With  w = U and  w = U q the  above  equations  may  be  rewritten, 
o 5 o 

respectively 


fiK  ’ 

- u 


[ i_iR»££v-*.,  j ^ ^/'{y  ^ 

+ L t;  (,ir^'i4nr).  f a’c’X  -{  |■^■]  4: 

^ r / ) 
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J 


I 

I 


where 


^ and 


S = 


d 


In  a rigorous  sense  the  nature  of  the  solution  for  coupled  axial 
and  pitch  motion  Is  determined  by  the  roots  of  a quartlc  character- 
istic equation,  and  the  roots  may  be  obtained  In  a variety  of  ways. 

The  vertical  plane  stability,  however,  may  be  more  easily  exam- 
ined by  considering  a simplified  form  of  Equations  (Cl)  and  (C2).. 

Under  the  assumption  that  * 0 and  the  cable  length,  h^.  Is  very 

large  with  respect  to  the  body  radius,  R , Equations  (Cl)  and  (C2) 

o 

can  be  decoupled  and  the  characteristic  equation  for  pitch  motion 
reduces  to  a quadratic  of  the  form 

as^  -I-  l>s  -h  C -c 

where 

For  the  case  of  the  Moored  Body  scheme  A,  a and  b are  both  negative 
(See  Table  IV)  therefore  the  sign  of  c determines  the  pitch  stability 
characteristics.  Roots  with  negative  real  parts  are  obtained  when 
c<  0,  which  means  that 

I.,R,CoQo^.  (C6) 


i 4 
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For  positive,  the  moored  body  Is  stable  at  low  current  speeds 


i 

» 

i 

' 

i 


l: 


i 


I 


but  unstable  at  higher  current  speeds.  For  negative;  stability  Is 
obtained  at  all  current  speeds. 

In  Equation  (C4)  It  may  be  noted  that  the  restoring  moment  provided 
by  the  vertical  center  of  gravity  acts  In  the  same  sense  as  the  cable 
tension  when  * 0 ®nd  thus  provides  an  additional  stabilizing  effect  to 
the  pitch  motion  as  expected.  Although  Equations  (C3)  and  (C4)  cannot 
rigorously  be  decoupled  when  *q  0,  an  approximate  condition  for  verti- 
cal plane  stability  which  Includes  the  effect  of  vertical  center  of 
gravity  but  neglects  any  surge  coupling  Is, 


HORIZONTAL  PLANE 

The  coupled  linear  equations  for  lateral  and  yaw  motion  of  the 
moored  body  In  a steady  current  are: 

Lateral  Force: 

nuir  Uof  + ) = (/^^i  y.  r -f  L V; 

- K.  (C8) 


Yaw  Moment: 

Ij  i-  -t-  in  C u -f  U.  <■ ) 

+•  A' r ^ +■ 


- A;,  r f i A',  i 

% c'  v;  r; 


Th«  above  aquation*  nay  b*  rawrlttan  raapactlvaly 


[ L'  ( Vj -m.')  s'  + L U.  V;,'s  - ] a 

[ L ) s'  + l U„<5  J % 

+ [ L'  C/v'/'Iz  )s  f i'U.  )]<•  = 0 


- Hi- 

/T71  - -7^ 


T — /p/  ></•  " 

IT  ^ ^ ^ TT 


and  s la  the  Laplace  operator,  . 

dt 

Equations  (CIO)  and  (Cll)  form  a pair  of  linear  homogeneous  differ- 
ential equations  for  1^2  ^ roots  of  the  characteristic  equa- 

tion for  the  above  set  provide  the  basic  form  of  the  solution. 

The  characteristic  equation  Is  a cubic  of  the  form 


where : 


CL  = [( ) -i)) -fri'ii ) ( A/; -m'*; ) 

- ( ) - V/  C ^0 ) ] r 
, _ li' 


f/j.u,  v:  L* 
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If  all  of  the  roots  of  Equation  (C12)  have  negative  real  parts  the 
body  la  stable  In  the  horizontal  plane.  For  a cubic  equation  a neces- 
sary and  sufficient  condition  that  all  of  the  roots  have  negative  real 


parts  is  that: 

(a)  all  coefficients  have  the  same  sign. 

(b)  be  - ad  > 0 


(C13) 


1 


1 

I 

I 
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ID  DTNSRDC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

(2)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMBERED  AS  TM  SERIES 
REPORTS,  NOT  FOR  GENERAL  DISTRIBUTION. 


